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Abstract 

A power law cosmology is defined by the cosmological scale factor evolving as t". In 
this work, we put bounds on a by using the joint test of the SNe la data from Supernova 
Legacy Survey (SNLS) and Hi^z) data with curvature constant A; = 0, ±1. We observe 
that the combined analysis with SNLS and H(z) data favours the open power law 
cosmology with a. = l.Slto^os- It is also interesting to note that an Einstein - de Sitter 
model (a = 2/3) is ruled out at 2a level. 

1 Introduction 

Our universe is very well explained by the Standard Cosmological Model (SCM) based on the 
Hot Big-Bang theory and the inflationary scenario. However, there are still some features 
of the universe which cannot be understood within the SCM. One of the major unsolved 
problem is the cosmological constant problem. The standard model fails to explain why 
the energy density of the vacuum is 120 orders of magnitude smaller than its value at 
the Planck time The problems in the SCM and the availability of precise data from 
various observations have encouraged cosmologists to explain the observed universe through 
alternative cosmological models. 
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One of the interesting alternatives is a Power Law Cosmology. In such a model there is 
a power law evolution of the cosmological scale factor, a{t) oc t". The power law evolution 
with a > 1 has been discussed at length in a series of earlier articles [2l[3l|4]. The motivation 
for such a scenario comes from the fact that it does not encounter flatness and the horizon 
problem at all. Another interesting feature of these models is that they easily accommodate 
high redshift objects and hence alleviate the age problem. These models are also purged of 
the fine tuning problem [5[ |6]. Such a scaling is a generic feature in a class of models that 
attempt to dynamically solve the cosmological constant problem [3|, HI O [6l [7] . 

A power law evolution of the cosmological scale factor with a ~ 1 is surprisingly an excellent 
fit to a host of cosmological observations. Any model that can support such a coasting 
presents itself as a falsifiable model as far as classical cosmological tests are concerned as it 
exhibits distinguishable and verifiable features. An evolution of this nature is supported by 
classical cosmological tests such as the galaxy number counts as a function of redshift and the 
data on angular diameter distance as a function of redshift [2]. However, as these tests are 
marred by evolutionary effects (e.g. mergers), they have fallen into disfavour as reliable tests 
of a viable model. With the discovery of Supernovae type la, SNe la, as reliable standard 
candles, the status of the Hubble test has been elevated to that of a precision measurement. 
The Hubble plot relates the magnitude of a standard candle to its redshift in an expanding 
FRW universe. We demonstrated that linear coasting cosmology accommodates the high 
redshift objects while the standard model could not [8l[9]. Such a model is also comfortably 
consistent with the gravitational lensing statistics [8] and the primordial nucleosynthesis [lOJ. 

The plan of the paper is as follows. In Section 2, we give the basic equations for open, 
closed and fiat power law scenarios. In Section 3, we find the constraints on the cosmological 
parameter a from a joint test of the Supernova Legacy Survey SNe la data set (SNLS) and 
the H(z) data. The joint test is performed for open, closed and fiat power law cosmologies. 
The results are summarized in Section 4. 

2 Power Law Cosmology 

For a FRW metric, the line element is 



Here k = ±1,0, is the curvature constant, t is the cosmic proper time and a{t) is the 
cosmological scale factor. 

The expansion rate of the universe is described by a Hubble parameter, H{t) = a/ a. The 
present expansion rate of the universe is defined by the Hubble constant Hq. Here and 
subsequently the subscript on a parameter refers to its present value. 



rfs^ = c^dt^ - a^{t) 



+ r^{de^ + sm^ed(l)^) 




1 — kr'^ 



2 



In this paper, we study a general power law cosmology with the scale factor given in terms 
of a dimensionless parameter a 

In this model, H{t) = a/t and Hq — a /to- The scale factor and the redshift at time t are 
related to their present values by 

1 + z . (3) 



a(z) t 

The present 'radius' of the universe is defined as 



ao = ^ . (4) 



The age of the universe at redshift z is given as 



The dimensionless Hubble parameter is defined as: 

E{z)^^ = {l + zy/-. (6) 

For the power law cosmology, the luminosity distance between two redshifts Zi and Z2 is 

<i,(,.,,,) = £(l±£5)5(_^{(i + ,,)^-(i + ,0=J^}) . m 



In the limiting case, a —>■ 1, we obtain 



ddz^, Z2) = [/n(l + Z2) - ln{l + z,)] . (8) 

-no 



Here 



S{x) — sinh(a;) for A; = — 1 
= sin(a;) for k — +1 
— X foT k — 

(9) 
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3 Observational Tests 



3.1 Constraints from the Supernova Legacy Survey SNe la data 
set (SNLS) 

Type la supernova (SNe la) are excellent cosmological standard candles for estimating the 
apparent magnitude m{z) at peak brightness after accounting for various corrections. In this 
work we use the SNLS data set of 115 SNe la data points with redshift z < 1 [12\ . 

For a standard candle of absolute magnitude M, the apparent magnitude m{z) can be ex- 
pressed as: 



m{z) = M + 5 logio Diiz) . 
Here Dl{z) is related to the luminosity distance: 

Ho 



Dl{z) 



-dUO^z) 
c 

:i + z)s 



a 



a — 1 



l + z) 



(10) 



and 



M = M-51ogio/i + 42.38 



(12) 



is the "zero point" magnitude. We use Hq = /ilOOKms ^ Mpc ^. 
The distance modulus, fi{z), is defined as 

ix{z) = m{z) - M = 5 logio Dl{z) - 5 log^o h + 42.38. (13) 



The chi-square function is defined as 

115 
i=l 

where fiexp is the expected distance modulus for a supernova at a given redshift z and cTj is the 
error due to intrinsic dispersion of SNe la absolute magnitude and observational uncertainties 
in SNe la peak luminosity. These errors are assumed to be Gaussian and uncorrelated. The 
observed distance modulus, fiobs, is given by the supernovae data set. 



fil (h,a,Zi) - filf^^iziY 



(14) 
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3.2 Constraints from H(z) data 



Simon, Verde and Jimenez (2005) used differential ages of passively evolving galaxies to 
determine the Hubble parameter as a function of redshift, H[z) [I3]. They use a sample 
of absolute ages of 32 galaxies taken from the Gemini Deep Deep Survey (GDDS) and the 
archival data to obtain 9 data points of H{z) with 0.09 < z < 1.75. The Hubble parameter 
and the differential age of the universe, dz/dt, are linked by the equation: 

The details of the method for calculating the dz/dt from the absolute age is given by Simon, 
Verde and Jimenez (2005) [13] . 

The H{z) for power law cosmology is given by: 

E{z) - ^ = (1 + ^)^/^ 

-no 

In order to put bounds on the model parameter, a, we define the quantity: 

X{h,a)=}_^ 2 (16) 

i=l 



Where H^xp is the expected value of the Hubble constant in the power law cosmology. Hobs 
is the observed value and is the corresponding la uncertainty in the measurement. 



3.3 Joint Test: SNe la + H(z) 

We find the constraints on the cosmological parameter a from the joint test of SNe la and 
H(z) data sets. In this joint test we define the quantity 

2 _ 2 I 2 n 7^ 

Xjoint — XsNe + Xh(z) ' \^ ' ) 

where x%Ne is given by Eq. (fT^ and x\{{z) by Eq. (fT6l) . 

Considering /i to be a nuisance parameter, we marginalize over h to obtain the probability 
distribution function defined as: 

L{a) = J e-^j°i-('^' P{h) dh. 
Here P{h) is the prior probability function for h which we assume to be Gaussian: 

= -f^- exp[-- -2 ], 

V27ra/, 2 Oy^ 
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Closed Flat 



Open 



Set A Prior 



Xl = 1-51 



'X.p 



1.28 



1.15 



2.: 



-0.23 ^ 
-0.19 (-^min 



1 «o+0.10 
1.62_o,09 ftr 



1.31 



h0.06 
-0.05 



Set B Prior 



xl = 1-53 



xl = 1-30 



xl = i.ir 



2.28^Q CHrnin 



1 62+°-^° rT ■ 



1.31 



hO.06 
-0.05 



Table 1: The best fit value of a and xl = Xmm/(degree of freedom) for the three models. 



where hots is the value of h (and is the error in it) as suggested by independent observa- 
tions. In this paper, we also study the effect of different priors on the result. We use two set 
of priors: 

(i) Set A: hobs = 0.68 ± 0.04 as obtained from the median statistics analysis of 461 mea- 
surements of Hq [14J. 



fii) Set B: h 



obs 



0.77 ±0.04 as suggested by the Chandra X - ray Observatory results 



The best fit model parameter is obtained by minimizing the modified (obtained after 
marginalization over h) : 



X 



-21nL(a) 



'1^ 



We performed the joint analysis for open, closed and flat power law cosmologies. The best 
fit value of a and the constraints on it seem to be independent of the choice of prior in all 
the three models. The result obtained with both the priors are summarised in Table 1. The 
joint analysis also shows that the best fit scenario is an open model with ami-^ 



1 O1+0.06 
-'-•'J-L-0.05- 



4 Summary 

In this paper, we study observational constraints on the power law cosmology, a{t) oc t". This 
model of the universe has very interesting features which makes it unique when compared 
to the other models of the universe. Firstly, for a > 1 it does not encounter the horizon. 
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flatness and age problem [21 [Sj H] . Secondly, such an evolution is a characteristic feature of 
models that dynamically solve the cosmological constant problem. Statistically this model 
may be preferred over other models as we have to fit only one parameter, a. 

In the work presented here, we use the joint test, which uses the SNLS data and H(z) data, 
to put constraints on the parameter a. To begin with, we work with all the three models 
- closed, fiat and open and put constraints on a in these three cosmologies. The results on 
the cosmological test are summarized in Table 1. Fig. 1-3 show variation of with a for 
the three models ( with set B prior). We also mark the parametric space allowed at 90% CL 
in the figures. We make the following observations: 

• For the three models under consideration, the value of x^, the best fit value of a and the 
constraints on it seem to be independent of the choice of prior. 

• The joint test favours an open power law cosmology with = l-Sllaos- As can be seen 
in Table 1, is minimum for an open power law cosmology. 

• The joint analysis does not rule out fiat and closed power law cosmologies. However, we 
do observe that the constraints on a are tighter for an open model. 

• The joint analysis rules out linear coasting [a = 1) in all the three cosmologies even at 90 
% CL. (see Fig.l, Fig. 2 and Fig. 3). 

These observations match the conclusions of Zhu et al. (2007) |16j. They test the power-law 
cosmology against the recent measurements of the X-ray gas mass fractions in clusters of 
galaxies. They conclude that the best fit is an open model (with amin = 1.14 ±0.05) though 
the fiat and closed models can not be ruled out. 

In the past, various observational tests have been used to put constraints on the parameter 
a in an open power law cosmology, such as gravitational lensing. Old High Redshift Galaxies 
(OHRG), SNe la and X-ray gas mass fractions in galaxy clusters. Constraints obtained from 
the other tests along with the constraints obtained from the SNe la data and H(z) data are 
summarized in Table 2. The interest in an open power law cosmology is on account of the 
fact that a whole class of dilaton gravity models that dynamically solve the cosmological 
constant rely on a vanishingly small effective gravitational constant in the early universe 
[71 [T7] . This gives an open FRW model for any reasonable equation of state of matter. 

We observe that the joint analysis of SNLS and H(z) data favours an open power law cos- 
mology with a > 1. We further observe that the joint analysis (SNLS + H(z)) done in this 
letter rules out the Einstein-de Sitter universe (a = 2/3). 

Since the joint test of SNLS and H(z) data presented in this work favours an open power 
law cosmology, for the sake of completeness we find bounds on a in open model separately 
using the SNLS data and the H(z) data. We once again marginalize over h to find xt 
the best fit values using each test. We find that: 

1. Constraints from SNLS data: For set A prior, we get the best fit value a = 1.4211q;q7 
with xl = 1-07 With set B, we get the same constraints on the parameter a as obtained 
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from set A but with xl = 1-09. We find that the constraints on a do not depend upon the 
choice of the prior. We, therefore, conclude that the SNLS data favours a > 1 (best fit value 
being a = 1.42lo;o7 ). This observational data rules out linear coasting cosmology {a = 1) 
even at at 2a level. 

2. Constraints from H(z) data: We find that the H(z) data provides tight constraints 
on the model parameter a. With set A prior we obtain best fit value a = 1.02]!^o;gg with 
xl = 0.834. With set B prior we get the best fit value as a = LOrl^og with xl = 1-06. We 
observe that for this test the constraints on a weakly depend upon the choice of priors. The 
H(z) data, however, strongly favours linear coasting cosmology with the best fit value. 

We summarize: An open power law cosmology with a > 1 is in excellent agreement with the 
present day observations. This makes it an attractive alternative. In fact, the possibility of 
an open linear coasting model as a viable model cannot be ruled out (as suggested by the 
H(z) data). Concordance of the power-law cosmology with CMB anisotropy measurement 
is a major area to be explored. There are large numbers of surveys that are ongoing or 
have been proposed. With the flood of new data (and the possibility that the observational 
techniques will be improved), the task ahead is to find models of the universe that can 
explain these observations. It will be interesting to investigate how the future observations 
will change the constraints on a. 
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Closed Power Law Cosmology 
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Figure 1: Results with Gaussian prior hobs = 0.77 ± 0.04 in a closed power law cosmology. 
The vertical lines at a = 1.98 and at a = 2.69 mark the parametric space allowed at 90% 
CL. The minimum of occurs at a = 2.28 . 
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Flat Power Law Cosmology 
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Figure 2: Results with Gaussian prior hobs — 0.77 ± 0.04 in a flat model. The vertical hues 
at a = 1.48 and at a = 1.79 mark the parametric space allowed at 90% CL. The best flt 
value occurs at a = 1.62. 
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Figure 3: Results with Gaussian prior hobs — 0.77 ±0.04 in an open model. The vertical hues 
at a = 1.23 and at a; = 1.41 mark the parametric space allowed at 90% CL. The minimum 
occurs at a = 1.31 . 
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